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ABSTRACT 

Limit loads under plane stress and plane strain are found for a 
circumferentially reinforced elastic-plastic ring subjected to interior 
pressure. These are used as bounds on an estimate of the failure pressure of 
a SiC/Ti test ring that is being fabricated and tested under the 
co-sponsorship of NASA Lewis Research Center and Pratt and Whitney 
Aircraft. The ring is to serve as a benchmark against which deformation and 
failure analysis methods can be assessed. An anisotropic perfect plasticity 
idealization of the SiC/Ti ring material is made and used in the limit load 
calculations. An estimate of the failure pressure of the NASA/PW 
benchmark test ring is given. 
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INTRODUCTION 


A plastic limit analysis is made for estimating the failure pressure of a 
SiC/Ti test ring subjected to a monotonically increasing pressure at 427 C. 
The circumferentially reinforced SiC/Ti ring models a critical aircraft 
engine component. The fabrication and testing of the ring is co-sponsored 
by NASA Lewis Research Center and Pratt and Whitney Aircraft as a 
benchmark for assessing the applicability and accuracy of structural 
analysis and failure prediction methods for metallic composites. 

Following Robinson and Duffy (1990) and Binienda and Robinson 
(1990), the SiC/Ti material is considered a pseudohomogeneous, locally 
transversely isotropoic continuum and is idealized as elastic perfectly 
plastic. The concepts of limit analysis of elastic-plastic continua are used 
under the assertion that strains large enough to cause rupture occur only as 
the plastic limit load P 1 is reached (fig. 1). This approach circumvents 
details such as residual stresses that affect the initial yield pressure Py 
(fig.l) and the intermediate elastic-plastic response but not the limit 
pressure P 1 . The failure pressure of the NASA/PW ring is bounded by 
calculating the limit pressures of the two limiting cases of plane stress and 
plane strain. 

The objective is to provide an estimate of the failure pressure of the 
NASA/PW benchmark ring using approximate methods useful to the 
design engineer without conducting a detailed and costly structural analysis 
on a large computer. Limit analysis provides closed form expressions for the 
failure pressure, or bounds on it, valid for a range of ring geometries and 
material properties; such simple and flexible tools are useful to the designer, 
especially in the early stages of design. The limit load calculations also 
provide a check for verifying the correctness of a more comprehensive finite 
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element analysis of the NASA/PW ring. 

The authors anticipate a sequel to this publication comparing the 
results obtained here with experimental results on the NASA/PW 
benchmark ring as they become available. 

LIMIT ANALYSIS 


We first state the anisotropic plasticity theory used in the analysis. 
Following Robinson and Duffy (1990), dj(xj) is a unit vector field 
designating the local fiber direction (normal to the plane of transverse 
isotropy) at each material point. Dij = djdj is a symmetric second rank 
orientation tensor. The yield function <J>(cij) is 

tfaij) = J2--|~(I -K , (1) 


written in terms of two invariants of the deviatoric stress Sij 

J2 = -5-SijSji 
and I = DjjSji 


and Dij, viz, 

( 2 ) 

( 3 ) 


The parameter 0 < ( < 1 designates the degree of anisotropy. ( = 0 
corresponds to full isotropy whereas ( =1 indicates unlimited strength in 
the fiber direction. K is the yield stress under transverse shear. 

The flow law is 

«ij = ATij (A > 0) if <J) = 0 and = 0 (4) 

cij = 0 if (j) < 0 or (J> = 0 and r^j < 0 
in which Ty } = sjj - -j-fl(3Dij - fy) (5) 
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Note that with ( = 0 this theory reduces to the classical J 2 (v.Mises) 
theory of a perfectly plastic solid. 

Reduction of eqs.(l)-(5) to uniaxial stress gives 

Y, = K V~37(HJ Y t = K/12/(W) (6) 


as the respective yield stresses along and transverse to the fiber direction. 

In terms of cylindrical coordinates appropriate for the ring (fig.2) the 
nonzero stress components are <7 z ,<r c and a T . The inner and outer radii are 
denoted as a and b, respectively. Circumferential fiber reinforcement 
dictates that D 22 = D c = 1, otherwise Dij = 0. 

The deviatoric stress components are 

S Z = -y^(<7 z -<7 r )+(<7 z -<7 c )] 

Sc = -jii <7c-*z)+( ff c~&r)] (7) 

3 r = “J"{( c)+( r — °z)] 

and the components of r*j axe 

T z = s z + "" 2 _ C 3 c 

r c = (1-C) Sc (8) 

r r = s r + “ 2 “(s c 

In these terms, the yield function eq.(l) becomes 

i| = y(>- , U ! )-T(3(-2K , -K 1 (9) 


and the flow law is 
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( 10 ) 


€ p = AT p if (j) = 0 and r p ij p = 0 

e p = 0 if ({) < 0 or <J> = 0 and TpO^ < 0 
with p = z,c,r 


The equilibrium equation relating <r c and <t t is 
d 

r - <V"0> 


and the compatibility equation in rate form is 



(11) 

(12) 


Generally, compatibility must be satisfied by the total strain rates. 
However, at the limit load where the elastic strain rates vanish, eq.(12) 
must be satisfied by the plastic strain rate field. 

We shall first consider the plane strain problem, i.e., e z = 0. Using 
this condition in eqs.(7),(8) and (10) results in 



Thus we can eliminate cr z in eqs.(7) and in the yield function eq.(9) giving 


or 


<j> = -^(<rc-<r r ) 2 - K2 = 0 


<Tc-0Y - 



K 


(14) 

(15) 


Fig.3 shows the yield surface, eqs.(15), in the dimensionless stress space 
(<r c /2K, <r r /2K). Also shown is a typical stress profile a-b for pressure 
loading under plane strain; evidently, the positive root in eqs.(15) is the 
physically meaningful one. Combining this and the equilibrium eq.(ll), 
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yields a separable differential equation in <r T and r. This equation can be 
integrated from r = a to r = b, making use of the boundary conditions 
<r r (a) = -P 1 and <r r (b) = 0, to yield a lower bound on the (dimensionless) 
limit pressure P 72 K for plane strain, i.e., 


P 1 _ 1 rs 

j H 




So far, only equilibrium, the constitutive law and the boundary conditions 
are shown to be satisfied ( not compatibility), restricting eq.(16) to be a 
lower bound on the limit pressure. 

The dimensionless equilibrium stress field at the limit pressure is 

(17) 

7 ) ( 18 ) 

with a z j 2 K given by eq.(13). Inserting this stress field in eqs.(7),(8) and 
( 10 ) results in the associated plastic strain rate field 




where, consistent with the perfect plasticity constitutive law, \ is taken as 
A(r) = 75 and A 1 = 2K ^ 


= const. 


This strain rate field is seen to satisfy compatibility, eq.( 12 ), thus elevating 
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eq.(16) to the exact limit pressure for plane strain under the present 
material idealization. 

We now consider the plane stress problem in which <r z = 0 . With this 
condition in eqs.(7), the yield condition eq.(l) becomes 

♦ = J KKk 2 - (KVr + (1- - K 2 (20) 

At yield, ({> = 0, we can solve eq.(20) for <r c giving 



The yield surface, eq.(20), is shown in the dimensionless stress space 
(<t c /2K ( 0 >/ 2 K) of fig.4 for ( = 0 (isotropy), ( = 0.977 ( used subsequently 
for the NASA/PW ring) and ( = 1 (the limiting case of infinite strength 
in the fiber direction). Also indicated in fig.4 is a typical stress profile a-b 
for a plane stress solution under pressure loading with ( = .977. 

As for plane strain, the positive root in eqs.(21) is the physically 
meaningful one. It can be used to eliminate <r c in the equilibrium eq. (11), 
again resulting in a differential equation in <r r and r only. Integrating, as 
earlier, over r = a to r = b and invoking the boundary conditions <r r (a) 
= -P 1 and <r r (b) = 0, we obtain 





as a lower bound on the (dimenionless) limit pressure P!/2K for plane 
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stress. 

The dimensionless equilibrium stress field for plane stress is 
accordingly 

rO 


J 0" r 

2K 

with (7 c / 2K given by eq. (21) and u z = 0. 

Here, it is not such a simple matter to find an associated plastic strain 
rate field and show compatibility. However, we recognize that the strain 
rate field, eqs.(19), although not associated with the equilibrium plane 
stress field (eqs.(23) and (21)) through the constitutive law, is nevertheless 
kinematically admissible under plane stress and thus serves for establishing 
an upper bound on the limit pressure. This upper bound is obviously the 
limit load corresponding to plane strain, i.e., eq.(16). 

Another upper bound on the limit load for plane stress is that 
corresponding to a local ( out of plane) failure at the inner radius r = a 
directly under the applied pressure. This upper bound is P!/2K = 1 , as is 
easily calculated using the kinematically admissible slip line field of fig.5. 
For this local mode to occur, the point on the stress profile corresponding to 
r = a must be at C on the yield surface (fig.4). Point C (<r c /2K = -1/2, 
<r r /2K = -1) is independent of ( and has an associated strain rate vector 

with components e c = 0 and e r = - e z . 

Fig.6 is a plot of Py2K vs. ( for a pressurized ring in plane stress 
with b/a = 1.8, showing the lower bound, eq.(22), and the upper bounds, 
eq.(16) and P72K = 1. For an isotropic material ( = 0, the limit pressure 
(fig. 6) is P1/2K 8 0.58. As the material is strengthened circumferentially, 


2 dq 

rm=? i-n 


= In- 
r 




>-l 


(23) 
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the limit pressure increases (P!/2K s 0 84 at ( = 0.6). The limit pressure 
continues to increase with increasing ( until a critical value is reached 
(here, ( a 0.8 ) and the failure mode becomes one of local ( out of plane) 
flow at r = a. Further circumferential strengthening is ineffective in 
increasing the limit pressure. As ( increases beyond the critical value, the 
failure mode remains one of local flow with P1/2K = 1; we refer to this 
condition as over -reinforcement. Analogous behavior is found by Lance and 
Robinson (1972) in other reinforced structures, e.g., beams, plates, etc. We 
note that the upper bound in fig.6 corresponding to eq.(16), i.e.,the exact 
limit load under plane strain, continues to increase with ( , becoming 
infinite as ( -> 1. Thus, over-reinforcement in the sense discussed here does 
not occur under plane strain conditions. 

The critical ( in plane stress for a given b/a is obtained from 
eq.(22) by setting Py2K = 1. This is illustrated in the plot of b/a vs. ( in 
fig.7. Alternatively, this graph can be interpreted as giving the critical b/a 
ratio at which the local (out of plane) failure mode takes effect for a given 
(. For thick rings with b/a > the critical ratio, the ring fails locally at 
r = a and P1/2K = 1. 

Fig.8 shows the yield surface in plane stress, eq. (20), for ( = 0.8. 
The stress profile a-b in fig.8 (solid curve) corresponds to the plane stress 
solution discussed above for b/a = 1.8. Note that point a (i.e., r = a) on 
the stress profile coincides with point C (cf.,fig.4), which, as remarked 
earlier, is a necessary condition for the local (out of plane) failure mode. The 
dotted curve in fig.8 represents a stress profile typical for b/a < 1.8 , i.e., 
for a thinner ring. 

We shall apply the results of these limit pressure calculations to the 
NASA/PW test ring in the following sections. However, first we 
characterize the SiC/Ti material using available experimental data. 
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CHARACTERIZATION OF SiC/Ti 


The dotted curves in fig. 9 represent longitudinal and transverse 
tensile data on SiC/Ti at 427 C. The solid curves are the perfect plasticity 
idealizations of these data used in the subsequent limit pressure 
calculations. The yield stress in the fiber direction Y| is taken as being 
coincident with the fracture strength in that direction, viz., Y| = 200 ksi. 
Two values are taken for the transverse yield strength, Yt = 30 ksi and Yt 
= 35 ksi delimiting a range of possible idealizations. 

The material parameters ( and K corresponding to these choices of 
yield stress are obtained through eqs.(6) and are shown in the following 
table (Y|, Yt and K are in ksi): 

Table 1 

Y, Y t ( K 

200 30 .983 15.04 

200 35 .977 17.56 

APPLICATION TO THE NASA/PW TEST RING 

The NASA/PW ring has a circumferentially reinforced SiC/Ti core 
covered with a layer of Ti cladding. Fig. 10 shows the dimensions of the 
cross section and the relevant radii. Limit load calculations are made for 
two ring geometries (a and b are in inches), viz., 



Table 2 


b 

a 

b/a 

3.21 

2.78 

1.155 

3.25 

2.74 

1.186 
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The first, b/a= 1.155 corresponds to the actual radii of the SiC/Ti core. 
This calculation ignores the presence of the cladding, taking the strength of 
the ring to result entirely from the reinforced core. The second calculation, 
b/a = 1.186 incorporates effective radii (b » 3.25 and a « 2.74) as a 
means of approximating the small contribution of strength of the Ti 
cladding. Here, the choice of radii amounts to enlarging the core dimensions 
by about one third of the cladding thickness. A more prudent choice must 
come from engineering experience based on experimented observation. 

Critical values of ( , delimiting local failure in plane stress (fig.7), are 
( = .992 and .989 for b/a = 1.155 and 1.186, respectively. As each of 
these values exceeds ( = .977 and .983, identified in table 1 for the 
NASA/PW ring, we are assured that local failure is not a consideration 
here. 

Plots of P72K vs. ( for b/a = 1.155 and 1.186 are shown in 
figs.(ll) and (12), respectively. In each case, the dotted (upper) curve 
represents the exact limit pressure for plane strain and the upper bound for 
plane stress (eq.(16)). As discussed earlier, the upper bound P^/2K = 1 
relating to local flow at r = a does not play a role here. The solid (lower) 
curve in each of the figs. (11) and (12) is the lower bound for plane stress 
(eq.(22)). Specific values of the limit pressure for plane strain P* and the 
lower bound for plane stress P* corresponding to ( and K in Table 1 
and the ratios b/a in Table 2 are calculated as (K, P^andPl areinksi): 


c 

K 

Table 3 
b/a 

P 1 

1 s 

PI 

.983 

15.04 

1.155 

2U 

28.8 

.977 

17.56 

1.155 

26.5 

29.0 

.983 

15.04 

1.186 

2&1 

34.2 

.977 

17.56 

1.186 

30.1 

34.4 
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Table 3 indicates bounds on the limit pressure for b/a = 1.155 as 
25.3 ksi < P* < 29.0 ksi, spanning ± 7% about an average of 27.15 ksi. 
Recalling that b/a = 1.155 implies the strength of the ring is attributable 
solely to the SiC/Ti core, we anticipate these bounds might tend to 
underestimate the strength of the NASA/PW ring. 

Bounds for b/a = 1.186 (in which limited strength of the cladding is 
accounted for through specifying effective radii) are 
28.1 ksi < P 1 < 34.4 ksi . These range ± 10% about an average of 
31.25 ksi . The lesser value 28.1 ksi corresponds to the lower bound for 
plane stress with a transverse yield strength of Yt = 30 ksi (fig.9); the 

greater value 34.4 ksi relates to plane strain with Yt = 35 ksi. The 

authors submit these bounds as their best estimate of the failure pressure of 
the NASA/PW benchmark ring at 427C. 

CONCLUSIONS 

The following accomplishments and conclusions are made in this 
study: 

• The limit load is found for a circumferentially reinforced 

elastic-plastic ring (cylinder) subjected to interior pressure and in a 

condition of plane strain. 

• Upper and lower bounds on the limit load are found for a reinforced 
ring (disc) under conditions of plane stress. 

• These closed form solutions apply for all values of the material 
parameters 0 < ( < 1 and K (designating the degree of reinforcement and 
the yield stress in transverse shear, respectively), and any b/a (ratio of 
outer to inner radius). The results provide simple and useful tools for the 
engineer designing fiber reinforced metallic rings. 
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• As ( -* 1 the limit pressure of the circumferentially reinforced ring 
in plane strain increases indefinitely. Conversely, the limit pressure of the 
ring in plane stress is bounded by its transverse (a matrix) strength, 
ultimately failing by local plastic ( out of plane ) flow directly under the 
applied pressure. 

• The limit load calculations are applied to the NASA/PW SiC/Ti 
test ring giving an estimated range of the failure pressure at 427 C of a 28 - 
34 ksi. 

• The calculated bounds on the limit pressure are useful in providing 
a check on more comprehensive finite element analyses. 
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RADIAL DISPLACEMENT 


Fig.l Schematic diagram of pressure vs. radial displacement. 
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Fig.2 Cylindrical geometry (z, c,r) of circumferentially reinforced ring, b 
outer radius, a = inner radius. 
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Fig.4 Yield surfaces for plane stress in <r r /2K,<r c /2K space. Surfaces are 
shown for ( = 0, ( = .977 and ( = 1. Typical stress profile a-b for 
pressure loading with ( = .977. 


17 




PM = 4[(i/v5)(A/2)(h)(K)] 
P'/2K = 1 


Fig. 5 


Kinematically admissible slip line field for local (out of plane ) flow 
at r = a. Upper bound calculation giving P l /2K = 1. 
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Fig.6 P72K vs. ( for a pressurized ring in plane stress with b/a = 1.8. 
Lower bound, eq.(22), and upper bounds, eq.( 16 ) and P!/2K = 1. 
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